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Abstract - Self-propelled particles in an external gravitational field have been shown to dis- 
play both an increased sedimentation length and polar order even without particle interactions. 
Here, we investigate self-propelled particles which additionally are bottom-heavy, that is they feel 
a torque aligning them to swim against the gravitational field. For bottom-heavy particles the 
gravitational field has the two opposite effects of i) sedimentation and ii) upward alignment of the 
particles' swimming direction. We perform a multipole expansion of the one-particle distribution 
with respect to orientation and derive expressions for sedimentation length and mean particle ori- 
entation which we check against Brownian Dynamics simulations. For large strength of gravity or 
small particle speeds and aligning torque, we observe sedimentation with increased sedimentation 
length compared with passive colloids but also active colloids without bottom-heaviness. Increas- 
ing, for example, swimming speed the sedimentation profile is inverted and the particles swim 
towards the top wall of the enclosing box. We find maximal orientational order at intermediate 
swimming speeds for both cases of particles with bottom-heaviness and those without. Ordering 
unsurprisingly is increased for the bottom-heavy particles, but this difference disappears at higher 
levels of activity and for very high activities ordering goes to zero in both cases. 



Introduction. — Actively moving particles, being in- 
trinsically out of equilibrium, exhibit a number of novel 
phenomena including interesting collective behaviour. Lo- 
cal alignment leads to the formation of moving swarms 
[IJ[2] or even turbulence [3] and hydrodynamic interac- 
tions may induce long-range order [4H6]. When particles 
slow down in the vicinity of others, this interplay between 
local particle density and speed leads to a motility induced 
phase separation [THT2"]. Autophoretic effects, where col- 
loids move along (self-generated) gradients of fields such 
as concentration or temperature, give rise to collective be- 
haviour reminiscent of chemotaxis in biological swimmers 

rani. 



tional order can be found even without particle interac- 
tions ghhh. 



But even in dilute suspensions, where particle interac- 
tions can be neglected, a number of interesting properties 
have been observed [T5HI5] . Activity changes the diffu- 
sive behaviour |19H22| and thus, for example, results in 
increased sedimentation lengths under gravitation [T5HT7] . 
While this effect can be mapped onto passive particles at 
higher temperatures, activity also aligns propulsion direc- 
tions of particles against an external field and orienta- 



Here, we investigate particles, which in addition to their 
activity exhibit bottom-heaviness, that is a particle's cen- 
tre of mass is displaced from its geometric centre in such a 
way that the particle tends to "point upward" and thereby 
swims against an external gravitational field. For strong 
gravity, bottom-heaviness dominates particle orientation 
but in the colloidal regime, where thermal diffusion is not 
negligible, we observe a combination of activity-induced 
and bottom-heavy orientational order. Starting from the 
Smoluchowski equation for particle position and orienta- 
tion, we derive an effective sedimentation process where 
activity and bottom-heaviness both enter into an increased 
effective diffusion constant and a reduced effective gravi- 
tation strength. We compare the analytically derived po- 
sition density profile and orientational distribution func- 
tion with results obtained from direct simulations of the 
equivalent Langevin equations. Furthermore, we identify 
a region in parameter space where the sedimentation pro- 
file is inverted and particles gather at the top wall of the 
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enclosing box and determine density profiles at small to 
moderate Peclet numbers. This is the regime relevant for 
colloidal systems used, for example, in [T4l[l"6] . 



Model system. — For bottom-heavy particles the 
centre of mass does not coincide with the geometric centre 
but is displaced toward the rear part of the particle, where 
front and rear are defined with respect to the swimming 
direction of the particle. The displacement is measured by 
the distance ro. Gravity g acts on the centre of mass and 
the torque r = — mr op x g tends to align the swimming 
direction p along — g. 

The resulting overdamped Langcvin equations for parti- 
cle position r and particle orientation p (|p| = 1) therefore 
read 
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g + V2Dt (1) 
(g x p) x p + y/2D~ r C x p (2) 



in the Stratonovich interpretation. The deterministic drift 
terms UoP and mg/7 give the swimming velocity of the 
active particle and the downward velocity due to gravity, 
respectively. The deterministic angular drift velocity is 
caused by the gravitational torque. While mg/7 actually 
requires the buoyancy mass Am = (p C oiioid — Pwatcr)-R 3 
(density mismatch with surrounding fluid) , the parameter 
in mro/7 r is indeed mass m. Here we assume that either 
Am w m or redefine the parameter ro as row/Am. 

Both translational and rotational motion are subject 
to stochastic velocities 



'2D£ and y/2D r £ where £ and 
£ each contain three independent Gaussian random vari- 
ables with zero mean and unit variance. In this Letter 
we assume that all noise is of thermal origin. So, the 
translational (D) and rotational (D r ) diffusion coefficients 
are connected with the respective friction coefficients via 
the Einstein relations D = ksT/"/ and D r = ksT/^/ r . 
For spherical particles rotational and translational fric- 
tion are related as 7,- = 4/3 R 2 j which, as the only source 
of randomness is thermal noise in the fluid, translates to 
D r = 3D/(4R 2 ) for the diffusion coefficients. 

We define a coordinate system by setting g = — ge z 
and non-dimensionalise spatial and temporal variables by 
setting r' = r/R and t 1 = tD/R 2 and omitting the 
prime. This introduces the following dimensionless num- 
bers, which we shall use as control parameters: the active 
Pcclct number or simply activity, Pc = vqR/D, the grav- 
itational Peclet number a = mgR/(k^T) 1 which is the 
particle radius over sedimentation length of the passive 
system, and excentricity or bottom-heaviness r' — r /R. 

The equivalent Smoluchowski equation for the one- 
particle density p(r, p, t) depending on both position and 
orientation then is 

d t p = —V ■ Jt — TZ ■ J r 

= V 2 p + (ae z - Pep) • Vp 

+ ^(K 2 p + ar K- {e z xpp)) (3) 



with 1Z = p x V p and the translational and rotational 
currents J t = (Pep — ae z )p — Vp and J r = — 3/4aro(e z x 
p)p — 3/ATZp. At the system boundaries the translational 
current has to vanish. 

Multipole expansion. — In order to derive an ap- 
proximative solution p(r, p,i) to Eq. in a systematic 
way, we use the multipole expansion 



p (u, (r,t) + p w {r,t) p 
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with coefficients p'°'(r,t), /oW(r,4), p' 2 '(r,£) depending 
on position and time but not on particle orientation. Since 
in our system of non-interacting scdimenting particles the 
only interesting dynamics takes place in z-direction, we 
assume that p only depends on z = r% and u = p% = cos 9. 
The multipole expansion for p reduces to an expansion in 
Legendre polynomials, 



p(z,u,t) = p {0) (z,t)+ p (L) (z,t)u 



+p [2) {z,t) ( u 
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(5) 



and the Smoluchowski equation ([3]) becomes 

d t p = d 2 p + (a - Pew) d z p 

+ ^ (d u (l - u 2 )d u p + ar a d u (l - u 2 )p) . (6) 

The position-dependent density n(z,t) = 
^_ 1 p{z,u 1 )du 1 which is an easy observable in exper- 
iments, is related to the zeroth-order expansion coefficient 
p {0) (z,t) simply as Q(z,t) = 2p(°\z,t). We begin with 
deriving an effective drift-diffusion process for density 
n(z,t) valid on time scales larger than the orientational 
dccorrclation time cx D^ 1 . It is of the form 



~d z J(l: 



(7) 



where Jo = a c ff £1 — D c ffd z Q is the effective density current 
with the drift and diffusive part. 

We now derive equations for coefficients p(°\ p^ and 
p( 2 ) from Eq. (j6|) making use of the orthogonality of Leg- 
endre polynomials. We first integrate Eq. ([6]) over u to 
obtain 



d t p^=d 2 JV+d z (ap 



Multiplying Eq. © by u and then integrating over u yields 
the equation for p^ 



dtpM = d 2 JV+ad zP W-PedJ-pW+pW 

\15 



(9) 
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The same procedure for u 2 — 1/3 gives the next order: 



V 2) = d 2 z pW+ad zP W--Ped z {^pW+pW 



(10) 



So, each equation couples p'*' to the next higher order 



Eqs. ((9]) and (| 10[) describe massive modes with 
<9t/9 ( - I ' ) ex p^\ i = 1,2, which relax faster than the slow 
mode dtp^ cx d z p^ from Eq. ([8]). We therefore omit the 
time derivatives in Eqs. (|9]) and (fTO]) while keeping the time 
dependence in Eq. ((5]) (cf. Refs. [TT1[TB]). Furthermore, 
we want to keep Eq. ([7]) at the drift-diffusion level with at 
most second-order derivatives and therefore drop second- 
order derivatives in Eq. ^ and second- and first-order 
derivatives in Eq. (fTU)) . We are still left with an equation 
for p( 2 ) that couples it to the next higher order p^ . This 
coupling is due to the bottom-heaviness of particles as 
would be the case with any rotational drift term. To close 
the relation, we assume p( 3 > <C and arrive at 
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(11) 



We plug this result into Eq. © and obtain an equation in 
p^\ Taking its derivative with respect to z and neglecting 
second-order derivatives gives d z p^ cx d z p(°\ which we 
re-insert into the equation for p' 1 ' to obtain 
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with A = (a 2 r - 2(ar ) 2 Pe/15)/(3 + (ar ) 2 /5). Using 
this result in Eq. dSJ, finally gives 



d t p^=D eS d 2 z p^-a eS d z p^ 



equal to the particle radius R meaning that particles ac- 
cumulate at the bottom wall. Then, Pe oc vq and thereby 
D e ff have to be sufficiently large, so that active particles 
can leave this densely packed colloidal phase. 

Within the regime of our approximations, we see that 
D c ff increases quadratically with Pe in accordance with 
the well-known result for the diffusion of (non-bottom- 
heavy) active particles [TW2"0] . The aligning torque (ar ) 
counteracts this trend as it tends to rectify particle swim- 
ming directions and thereby decrease the randomness in 
the active motion. For very large ar® the diffusion coef- 
ficient reverts to the passive case of purely thermal noise 
as swimming becomes fully deterministic. The effective 
gravitational strength a e g goes roughly linearly in a. Both 
bottom- heaviness (ro) and activity (Pe) work against grav- 
ity by causing the particles to swim upward. For a e g < 0, 
activity together with bottom-heaviness overcomes gravity 
and particles swim towards the top resulting in inverted 
density profiles. For a e g = the two effects balance and 
the density profile becomes uniform or constant in z (far 
from the confining walls). The limits of our approxima- 
tion again become visible for very large ro where Eq. f|15[) 
predicts a e g to tend to a instead of a — Pe. 

The swimming speed necessary to cancel sedimentation 
and to obtain a uniform density profile can also be cal- 
culated independently of the multipole expansion. A so- 
lution of Eq. © that is constant in z corresponds to the 
Boltzmann distribution of the particle orientation under 
the aligning torque, which gives p{z 1 u) cx exp(arou), also 
known as Mises-Fisher distribution. Since the bounding 
walls enforce that the average translational drift current 
along the z direction is zero, one has = Pe(u) — a. Cal- 
culating (u) from the Boltzmann distribution, we obtain 



(13) 
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with effective diffusion coefficient 

2Pe /Pe 



D, 



ft 



1 



3 + (ar ) 2 /5 V 3 



and effective gravity 

a c ff = a I 1 



r Pe 



3+ (ar ) 2 /5 



(14) 



We have thus identified the coefficients for the drift- 
diffusion equation for Vt{z,t) = 2p ( ~°\z,t) in Eq. (0. 

The steady state solution for f2(z) is an exponen- 
tial profile with inverse sedimentation length a e e/D c g, 
fl oc exp(— a e ftz/D e ff). We are now in a position to 
judge when the approximation to neglect derivatives of 
p^\ which we made in the derivation, is justified. As 
cx exp{-a cS z/D c ff) for i = 0,1,2 (Eqs. (fTT| [12] ITS])), 
this approximation is possible when the absolute value of 
the reduced inverse sedimentation length |a e ff/-D e ff| < 1- 
For awl passive colloids have a sedimentation length 



for the activity Pe* which cancels gravity so that a uniform 
density profile occurs. The approximate result for Pe* 
found from Eq. (|15|) by setting a c s = amounts to 



Pe* 



3 



a 2 r 



(17) 



^5) It corresponds to the first two terms of the Laurent series 
of Eq. (|16[) in org. 



Sedimentation length and orientational order. 

We solved the Langevin equations ^ and ([2]) for the sys- 
tem of bottom-heavy active particles by performing Brow- 
nian Dynamics simulations. We chose a cubic simulation 
box of size L 3 with L = 60 (in units of particle radii) with 
hard walls at z — and z = L and periodic boundary con- 
ditions in x- and y-direction. For low levels of activity (see 
Fig. [1] (a) for Pe = 3) the sole effect of bottom-heaviness 
is an increase in sedimentation length, while at higher ac- 
tivity, Pe = 8, bottom-heavy particles overcome gravity 
and produce inverted sedimentation profiles with particles 
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Fig. 1: (Colour online) (a) Spatial density profile and 
(b) orientational distribution uj(u) (within the bulk). Numer- 
ical data (dots) are computed from Brownian Dynamics simu- 
lations of Eqs. |T]) and ([2j|. Gravitational strength a — 1 and 
bottom-heaviness ro = 0.5 are kept constant, activity is Pe = 3 
(blue dots), Pe = 6.1 ~ Pe* (red dots) and Pe = 8 (green dots). 



swimming towards the top of the simulation box. The con- 
stant profile is realised for Pe = Pe* ps 6.1. All three cases 
are captured very well by Q(z) oc exp(— a e ff/D e ff) as can 
be seen in the comparison of simulations (dots) and theory 
(dashed lines) in Fig. Q] (a). The only significant devia- 
tions between the two appear at the hard walls in z which 
have been neglected in the theory. In the simulations, the 
boundary conditions in z give rise to an accumulation of 
particles swimming against the wall and thus an increase 
in density next to the wall. The net orientation at the wall 
also differs from the bulk (not shown) . This effect has been 
discussed in [T7] and accumulation of colloids next to the 
walls has also been observed in experiments [16]. 

The stationary orientational distribution function 
u)(u\z) = p(z,u)/Q(z) is found to be independent of po- 
sition z except in the vicinity of the bounding walls. We 
therefore write lo(u\z) = ui(u) for short. The orientational 
distribution also depends on the particles' activity as can 
be seen from Fig. [1] (b). This ordering is due to two ef- 
fects: Firstly, the bottom-heaviness simply aligns particles 



against the vertical while the second effect due to activ- 
ity works more subtly and in conjunction with the spa- 
tial density profile as shown in Ref. [17] • An excess of 
particles at position z — Az will mean that more parti- 
cles reach position z swimming upwards from below than 
swimming downwards from above and particles at z will 
predominantly point upwards. At low activity, Pe = 3, the 
sedimentation profile is more pronounced and it produces 
alignment in the same direction as bottom-heaviness. In 
contrast, an inverted profile such as for Pe = 8 works 
in the opposite direction; it reduces orientational order 
and results in a flatter distribution, see Fig. [1] (b). The 
orientational distribution u(u) can also be derived from 
the multipole expansion Eq. ([5]). After insertion of p^ % \ 
i = 0, 1, 2 we find 



with 

ui(u) cx 1 + 



p(u, z) = fl(z)ui(u) 



1 



(18) 



. ^(Pc - 3 A) + ^ar 
3+ (m-o) 2 /5 VArfr 2 



2u + aro ( u 



1 



(19) 



As we only performed the expansion up to second order, 
ui{u) is quadratic in u by construction. A better fit of the 
data, but a less systematic result, can be found from a 
perturbation ansatz in a, ui(u) = exp(a/i(u) + a 2 /2(u)), 
where fi(u) and f2(u) can be determined from inserting 
Eq. (|18p into Eq. ((6]). Both results are plotted against 
simulation data in Fig. Q] (b). 

We find that sedimentation lengths increase with activ- 
ity Pe for both bottom-heavy and non-bottom-heavy par- 
ticles. However, the increase for bottom-heavy particles is 
stronger as Fig. [2] (a) demonstrates where we plot inverse 
sedimentation length a c g /-D e ff over Pe. For large Pe both 
curves converge to a c tf/D c g = (infinite sedimentation 
length) as the effective diffusion diverges quadratically in 
Pe, that is the contribution of activity to random mo- 
tion increases faster with Pe than the contribution to di- 
rected motion. However, the curve for bottom-heavy par- 
ticles does so from below as for large Pe upswimming wins 
over sedimentation and the drift term becomes negative, 
a e s < 0. This regime of inverted profiles is visible even 
better in Fig. [2](b), where inverse sedimentation length is 
plotted against gravitational strength a. For intermediate 
values of a, bottom-heavy particles are oriented upwards 
but gravity is not yet strong enough to pull particles down 
and we find maximally inverted profiles for a w 3 (dashed 
blue curve). For non-bottom- heavy particles (solid red 
curve) the inverse sedimentation length increases linearly 
for small a, in agreement with Ref. [17]. Deviations be- 
tween simulation and theory appear once sedimentation 
lengths become too small, as discussed before. Parameter 
ranges relevant for active colloidal systems in Ref. [TB] are 
a = 0.1 and Pe = to Pe = 6 but larger values could 
be realised by increasing particle radii since a oc i? 4 and 
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Fig. 2: (Colour online) Inverse sedimentation length a e ff/-D e s 
(a) over activity Pe at constant a = 1 and (b) over gravi- 
tational strength a at constant Pe = 8. Lines are theoreti- 
cal predictions, dots are simulation results. Curves for non- 
bottom-heavy particles (ro = 0) are shown in solid red, curves 
for bottom-heavy particles (ro = 0.5) are dashed blue. 



Fig. 3: (Colour online) Mean orientation (u) (a) over activ- 
ity Pe at constant a = 1 and (b) over gravitational strength 
a at constant Pe = 8. Lines are theoretical predictions, dots 
are simulation results. Curves for non-bottom-heavy particles 
(ro = 0) are shown in solid red, curves for bottom-heavy par- 
ticles (ro = 0.5) are dashed blue. 



Pe cx R. The region of inverted sedimentation for bottom- 
heavy particles is especially interesting as tuning of a in a 
centrifuge might be used to separate bottom-heavy from 
non-bottom-heavy colloids. 

As has been mentioned above, the orientational distri- 
bution changes with activity Pe and becomes flatter for 
large values. Indeed, we find a maximum in mean particle 
orientation (u) = ft(z)~ 1 J_ ± p(u, z)udu for intermediate 
values of activity. From (u) = p^ 1 ' /(3p( ') we then find 

(„) = _i_ (**(*1-a\+ ±ar ) (20) 
V ' 3 + (w ) 2 /5 \D eS \ 3 ) 2 7 V 1 

which at large Pe can be roughly viewed as an activity- 
induced part cx a e ff/D e s (although depending on the spa- 
tial profile which is changed by the presence of bottom- 
heaviness) and the direct contribution due to bottom- 
heaviness. In Fig. [3] (a) for fixed a = 1, polar order is max- 
imal for Pc ~ 2 for non-bottom-heavy particles and for 
slightly lower values for bottom-heavy particles. Unsur- 
prisingly, ordering is stronger for bottom-heavy particles 
as bottom-heaviness itself directly contributes to particle 
alignment. For very large activities, however, (u) tends to 
zero in both cases. For non-bottom-heavy particles this is 



simply explained by a e g / D c g going to zero as Pe -2 . Phys- 
ically, this means the spatial density profile Cl(z) becomes 
constant and activity- induced order cx Pe a e s/D e g van- 
ishes. For bottom-heavy particles a e g/D e g goes to zero 
only as Pe -1 and activity- induced order becomes constant 
and cancels the bottom-heavy contribution. 

Ordering increases linearly with a for small a, see 
Fig. El (b), but flattens out for increasing values. For 
bottom-heavy particles the effect from increasing a is two- 
fold as it enhances both the passive alignment in an exter- 
nal potential and activity-induced alignment by sharpen- 
ing the sedimentation profile. Mean orientation therefore 
increases with a faster for bottom-heavy particles than it 
does for non-bottom-heavy particles although at the rela- 
tively large activity of Pe = 8 the curves indicate that the 
main effect for orientational order comes from the sedi- 
mentation profile. 

Conclusion. — To conclude, we performed a multi- 
pole expansion up to second order in orientation for the 
one-particle density of active bottom-heavy particles and 
determined a solution of the Smoluchowski equation in di- 
lute suspensions which is prerequisite to the study of inter- 
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acting particles. This approximation is useful in the col- 
loidal regime where active motion and gravity are not too 
dominant compared with thermal diffusion and the regime 
explored here is experimentally accessible, e.g., by chem- 
ically powered Janus particles [TUDiS. Comparing ori- 
entational order in bottom-heavy and non-bottom-heavy 
particles, we elucidated the influence of both activity (cf. 
Ref. [17]) and of bottom- heaviness. 

We derived analytical expressions for the sedimentation 
length of bottom-heavy particles and discussed the two 
regimes of increased sedimentation lengths [15UT7] and of 
inverted sedimentation profiles. For fixed gravitational 
strength a or fixed activity Pe, we find optimal values 
for the other parameter to maximise inversion. It is this 
density inversion and its instability under density fluctua- 
tions that, in the case of interacting particles, gives rise to 
pattern formation by bioconvection [2"5H2"5] (albeit in a pa- 
rameter regime of large Peclet numbers where thermal dif- 
fusion is small compared with active motion). We are cur- 
rently working on how this instability is induced by hydro- 
dynamic interactions between individual swimmers. Re- 
cently, it was demonstrated that in case of two Volvox al- 
gae hydrodynamic coupling leads to intricate bound states 
[26] . Finally, we note that alignment of particles can also 
be achieved by mechanisms other than bottom-heaviness. 
For example, for superparamagnetic colloids [27] or mag- 
netic Janus spheres |28] alignment is easily tuned by an 
external magnetic field whereas for biological swimmers 
chemotaxis or phototaxis result in similar phenomena |29) . 
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